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Abstract 

We report a study of the suppressed B~ —> and favored B —> Dtt~tt + 'k~ 

decays, where the neutral D meson is detected through its decays to the K^ir^ and 
CP-even I\ + K~ and 7r + 7r - final states. The measurement is carried out using a 
proton-proton collision data sample collected by the LHCb experiment, corresponding 
to an integrated luminosity of 3.0fb _1 . We observe the first significant signals in 
the CP-even final states of the D meson for both the suppressed B~ —> DK~tt + tt~ 
and favored B~ — > Dtt~tt + tt~ modes, as well as in the doubly Cabibbo-suppressed 
D —> K + ir~ final state of the B~ — > decay. Evidence for the ADS 

suppressed decay B~ — > DK~tt + tt~, with D — > K + n~, is also presented. From 
the observed yields in the P” — > DK~tt + tt~ , B~ — > ~ and their charge 

conjugate decay modes, the most probably value of the weak phase 7 corresponds to 
7 = (74^^g)°. This is one of the most precise single-measurement determinations of 
7 to date. 
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1 Introduction 


The study of beauty and charm hadron decays provides a powerful probe to search for 
physics beyond the Standard Model that is complementary to direct searches for new, 
high-mass particles. In the Standard Model, the flavor-changing charged currents of quarks 
are described by the 3x3 unitary complex-valued Cabibbo-Kobayashi-Maskawa (CKM) 
mixing matrix 11 , 21 , whose elements, U)j (i = u,c,t and j = d,s,b), quantify the relative 
i to j coupling strength. Its nine matrix elements can be expressed in terms of four 
independent parameters, which need to be experimentally determined. 

In general, decay rates that involve the i to j quark transition are sensitive to the 
magnitudes of the CKM matrix elements, V t] \ . The (weak) phases between different 
CKM matrix elements can be probed by studying the interference between two (or more) 
decay amplitudes. Particle and antiparticle amplitudes are related by the CP operator, 
where C signifies charge conjugation, and P refers to the parity operator. Under the 
CP operation, weak phases flip sign, leading to different decay rates for particles and 
antiparticles, if the weak and (CP-invariant) strong phases differ between the contributing 
amplitudes. Precision measurements of the magnitudes and phases of the CKM elements 
provide constraints on many possible scenarios for physics beyond the Standard Model. 

One of the least well-measured phases is 7 = arg[— (V u dV* b )/(V c dV* b )], which can be 
probed by studying the interference between b to u and b —> c transitions. The most 
promising method to determine 7 is to study the interference between B~ —> D°K~ 
and B —>■ D°K~ decays, when states accessible to both the D° and D° mesons are 
selected. These modes are particularly attractive for the determination of 7 because 
their amplitudes are dominated by only a pair of tree-level processes, leading to a small 
theoretical uncertainty |3|. Hereafter, we use D without a charge designation when the 
charm meson can be either a h° or h°. A number of methods, depending on the D 
decay mode, have been discussed in the literature, and are often grouped into three 
categories: (i) CP eigenstates, such as D —> K + K~ and D —> 7 r + 7 r - decays |4j[5] (GLW); 
(ii) flavor-specific final states, such as the Cabibbo-favored (CF) and doubly Cabibbo 


suppressed (DCS) D —> K ± tt t decays 


(ADS); and (iii) multi-body self-conjugate 


final states, such as D —)■ Kgir + n~ [| 8 | (GGSZ)Q 

Beyond this simplest set of modes, these techniques_are also applicable to modes with 
vector mesons, such as B 
6 -baryon decays, e.g 


to D*K~, P° to DI< 


*0 


A° b to DA 


10 


12 


multi-body final states of the recoiling strange quark system could be useful 13 


and B ^ to Deft [5J, as well as 
decays. It has also been suggested that other 

due to 


the larger branching fractions to these final states, and potentially a larger interference 
contribution. 

The current experimental measurements, averaged over several decays modes, are 


7 = (73_i 0 )° by the LHCb collaboration 14 ,7 = (691 \l)° by the BaBar collaboration 

£ 12 )° 


and 7 = ( 68 ^( 4 )° by the Belle collaboration |1 
fit to direct measurements of 7 is about 7° 17 


15 


The overall precision on 7 from a global 
To improve the overall precision on 7 , it 


1 The letters in the brackets are commonly used to refer to these general approaches, after the original 
authors. 
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is important to study a wide range of final states. 

In this article, we present the first ADS and GLW analyses of the decay B~ —y DX~, 
where the D meson is observed through its decay to , K + K~ and 7 r + 7 r^ final states 

and X~ = K~TT + n^. When specific charges are indicated in a decay, charge conjugation 
is implicitly included, except in the definition of asymmetries discussed below. The 
measurements use proton-proton (pp) collision data collected by the LHCb experiment, 
corresponding to an integrated luminosity of 3.0fb _1 , of which l.Ofb -1 was recorded at a 
center-of-mass energy of 7TeV and 2.0 ft> at 8 TeV. 


2 Formalism 


The formalism that was developed to describe the B 
the B~ —> DX~ case with only minor modifications 
states can be expressed as 
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—* DK modes can be applied in 
The decay rates in the CP final 


r (B -7- [h h + ) D X s ) oc 1 + r| + 2 nr B cos(S B - 7 ), (1) 

r(£? + -7- [h~h + ] D X+) oc 1 + r| + 2kv b cos(S b + 7 )- (2) 

Here, /i =b = 7 r ± or K ± , and [h~h + } B indicates that the state in brackets is produced in 
the decay of the neutral D meson. The quantities r B = \ A(B~ —* [h~h + ]^ 0 X~) / A(B~ —* 
[h~h + ] D oX~)\ and S B are the amplitude ratio and strong phase difference between B~ — * 
D° X~ and B~ —> D°X~ contributions, averaged over the DX~ phase space. The 

parameter k is a coherence factor that accounts for a dilution of the interference due to 

the variation of the strong phase across the phase space; its value is bounded between 
0 and 1. In principle, k can be obtained in a 1110 del-dependent way by a full amplitude 
analysis of this decay. Here, we consider it as a free parameter to be determined in the 
global fit for 7 . The strong parameters, r B , d B and k are specific to this decay, and differ 
from those obtained from other B —y DK modes. 

The decay rates for the D —» Jl ± 7r =F final states can be written as 


r (B~ -7- [K + tt~] d X~) oc r| + r 2 D + 2 nr B r D cos (5 B + 5 D - 7 ), (3) 

T(B + -7- [K~n + ] D X+) oc r| + r 2 D + 2nr B r D cos (5 B + S D + 7 ), (4) 

r (B~ -7- [K~n + ] D X~) oc 1 + ( r B r D ) 2 + 2 nr B r D cos (5 B — S D ~'y), (5) 

r (B + ->• [K + n~] D X+ ) oc 1 + (r B r D ) 2 + 2 nr B r D cos(d B - S D + 7 ). ( 6 ) 


Here, additional parameters r B and 5 B enter, which quantify the ratio of the DCS to CF 
amplitude, A(D° —>■ K + tt~)/A(D° —> K~tt + ) = r B e l5D . Values of r B and So are taken 
from independent measurements [18 , 

The determination of the CP observables in the /V —» DX~ decay uses the favored 
B~ —y DTt~TT +r K~ decay for normalization, denoted here as B~ —>■ DX^. For brevity, we 
will use X~ to refer to either Xj or X~. In addition, D Kn is used when both charge 
combinations are considered. 
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The CP observables of interest for the GLW analysis are the charge-averaged yield 
ratios 

h+h - = T(B~ -> [h+h-] D X;) + T(B+ -» [h+h~] D X+) 
cp+ r (B- -)■ [A'-7T+] D x s -) + r(A+ [A'+7T-] D x+) 

= 1 + + 2nr B cos5b cosy. (7) 


Because of the different D final states in Eq. [TJ systematic uncertainty due to the precision 
of the D branching fractions and the different selections is incurred. Following Ref. 
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we neglect CP violation in the B —y DX, and the favored D final state of B —» DX~ 


decays, and approximate R^J]_ by the following double ratio 


Rcp+ — 


Dh + h 

R S /d 

tdK-k ’ 
K s/d 


( 8 ) 


where 

. +k - = r(B- -> + r(B+ -» {hPrhxi) 

•!* r (B- -» [h+h~] D Xj) + r(B+ -> [h+h-] D X+y 

„ K , _ r(B^ -> [K-x+] D Xr + T(B* -» [A-+ I r-] D JC+) 

r(B- [A'-ir+] D A7 + r(B+ -> [A'+ir-] D .Y+) ■ 


(9) 

( 10 ) 


This double ratio has the benefit that almost all systematic uncertainties cancel to first 
order. The neglected CP-violating contribution of magnitude nr B \V us V c d/V uc iV cs \ < 0.01 is 
included as a source of systematic uncertainty. 

We also make use of the charge asymmetries 


Af 


r (b- f D x~) - r(R+ ->• f D x+) 
r (b- f D x-) + r(R+ f D x+) 


2 nr B sin S B sin 7 / Rcp+ , 


( 11 ) 


where / refers to either K + K , 7r + 7r or the CF K 7 r + final state in the D meson decay. 
For simplicity, small contributions from direct CP violation in D —> 7r + 7r _ and D —> K + K~ 


are not included here, but are accounted for in the fit for 7 14 


For the ADS modes, we measure the relative widths of the DCS to CF decays, separated 
by charge, as 


R x± = 


T(B ± -y [K t tt ± ] d X ± ) _ r B + r B + 2 nr B r D cos (5 B + 5 D ± 7 ) 
r(R ± ->■ [K ± tc^] d X ± ) 1 + r 2 B r 2 D + 2 nr B r D cos(<5 S - ± 7 ) ' 


( 12 ) 


The nearly identical final states in these ratios lead to a cancellation of the most significant 
sources of systematic uncertainty. Corrections to R x for D {) — D° mixing 120 are omitted 
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for clarity, but are included in the fit for 7 

All of the above equations, except for Eqsr |8 - 10 can be applied to either B ± —y DXf or 


B ± —y DX^ decays. The values of r B , S B and k differ between the favored and suppressed 
decays: however 7 is common to both. Most of the sensitivity is expected to come from 
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d )/A(B- -> D°XJ) is 0{\ 2 ) 
), where A = 0.2253 ± 0.0014 


as compared 
is the sine 
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the B ± —» DXf decays, since A(B —> D°X, 
to 0(1) for A(B~ ->• D 0 X~)/A(B~ D°X 
of the Cabibbo angle. Taken together, the observables that contain the most significant 
information on 7 are Rcp+, •A t x s h an d R Xa ■ Measurements of these four quantities 
constrain re, 5b, k and 7 . 

The product branching fraction for B~ —y DX~ decays, with D —>■ h + h~, is at the level 
of about 10 -6 . The small branching fractions, combined with a total selection efficiency 
that is of order 0.1%, makes the detection and study of these modes challenging. The 
corresponding ADS DCS decay mode is expected to have a yield of at least 10 times less 
than the CP modes, and is very sensitive to the values of re, 5b, k, and 7 (see Eqs. [3] 
and [d]). For this reason, the signal region of the ADS suppressed decays (both B~ —y DX^ 
and B~ —y DX~) was not examined until all selection requirements were determined. 


3 The LHCb detector and simulation 


The LHCb detector 22 is a single-arm forward spectrometer covering the pseudorapidity 


range 2 < 7 < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high-precision tracking system consisting of a silicon-strip vertex detector sur¬ 
rounding the pp interaction region, a large-area silicon-strip detector located upstream of 
a dipole magnet with a bending power of about 4 Tin, and three stations of silicon-strip 
detectors and straw drift tubes 1231 placed downstream of the magnet. The combined 
tracking system provides a momentum measurement with a relative uncertainty that varies 
from 0.5% at low momentum, p. to 1.0% at 200GeV/c, and an impact parameter mea¬ 
surement with a resolution of about 20 pm 24 for charged particles with large transverse 


momentum, pp- The polarity of the dipole magnet is reversed periodically throughout 
data-taking to reduce asymmetries in the detection of charged particles. Different types of 
charged hadrons are distinguished using information from two ring-imaging Cherenkov 
detectors 25 . Photon, electron and hadron candidates are identified by a calorimeter sys¬ 


tem consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter 
and a hadronic calorimeter. Muons are identified by a system composed of alternating 
layers of iron and multiwire proportional chambers 26 . Details on the performance of the 
LHCb detector can be found in Ref. 


27 


The trigger 28 consists of a hardware stage, based on information from the calorimeter 


and muon systems, followed by a software stage, which applies a full event reconstruction. 
The software trigger requires a two-, three- or four-track secondary vertex with a large 
Pt sum of the tracks and a significant displacement from all primary pp interaction 
vertices (PVs). At least one particle should have pr > 1.7GeV/c and xjp with respect 
to any PV greater than 16, where xjp is defined as the difference in y 2 of a given PV 
reconstructed with and without the considered particle. A multivariate algorithm 
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is 


used for the identification of secondary vertices consistent with the decay of a 6 -hadron. 

Proton-proton collisions are simulated using PYTHIA [30 with a specific LHCb con¬ 
figuration |31|. Decays of hadronic particles are described by EvtGen 32 , in which 
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33 . The interaction of the generated par- 


final-state radiation is generated using PHOTOS 
tides with the detector, and its response, are implemented using the Geant 4 toolkit 
as described in Ref. 
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35 . In modeling the B~ 


—> DX decays, we include several resonant 
and nonresonant contributions to emulate the X~ and XJ systems, as well as contributions 


from orbitally excited D states, e.g _Di(2420)° —» D°n + 7r~. The contributions are set- 
based on known branching fractions, or tuned to reproduce resonant substructures seen in 
the data. 


4 Candidate selection 


Candidate B decays are reconstructed by combining a D —* 
D —» 7r + 7r - candidate with an X~ candidate. A kinematic fit 
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Kn, D -a K+K - or 
is performed, where 

several constraints are imposed: the reconstructed positions of the A" - and B~ decay 
vertices are required to be compatible with each other, the D candidate must point back to 
the B~ decay vertex, the B~ candidate must have a direction consistent with originating 
from a PV in the event, and the invariant mass of the D candidate must be consistent 
with the known D° mass 21 . The production point of each B~ candidate is designated 
to be the PV for which the Xip is smallest. 

Candidate D mesons are required to have invariant mass within 3<7 d (2.5ctd for 
D —>• 7r - 7r + decays) of the known value, where the mass resolution, crp, varies from 
7.0MeV/c 2 for D —> K + K~ to 10.2MeV/c 2 for D —>• 7r + 7r - decays. Unlike the D mesons, 
the invariant mass of the X~ system covers a broad range from about 0.9 — 3.3GeV/c 2 . 
Candidates are required to have an invariant mass, M(X ~) < 2.0GeV/c 2 . For the X~ 
system, we also require the K~n + invariant mass to be within 100MeV/c 2 of the known 
K *° mass. The latter two requirements not only improve the signal-to-background ratio, 
but should also increase the coherence factor k in the final state. 

To improve the signal-to-background ratio further, we select candidates based on particle 
identification (PID) information, and on the output of a boosted decision tree (BDT) 137, 


38 classifier. The latter discriminates signal from combinatorial background based on 
information derived primarily from the tracking system. For the BDT, signal efficiencies are 
obtained from large samples of simulated signal decays. Particle identification efficiencies 
are obtained from a large D* + —> D°tt + calibration data sample 25 , reweighted in pp, V 


and number of tracks in the event to match the distributions in data. The effect of the 
BDT and PID selection requirements on the background is assessed using sidebands well 
away from the B~ peak region. In the optimization, a wide range of selection requirements 
on the PID and BDT outputs are scanned, and we choose the value that optimizes the 
expected statistical precision of the B~ —> DX~ signal yield. Expected signal yields 
are evaluated based on known or estimated branching fractions and efficiencies obtained 
from simulation (for the BDT) or D* + —>■ D°n + , D° —>■ K~n + calibration data (for the 
PID). Due to the smaller expected yields in the ADS modes, separate optimizations are 
performed for the GLW and the ADS analyses. Using simulated decays, we find that the 
relative efficiencies for B~ —» DXj and B~ —>■ DXt decays across the phase space are 
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compatible for the GLW and ADS selections. Due to the uniformity of the selections, and 
the fact that the observables are either double ratios, e.g. Rcp+, or ratios involving almost 
identical final states, the systematic uncertainty on the relative efficiencies is negligible 
compared to the statistical uncertainty. 

Several other mode-specific requirements are imposed to suppress background from 
other 6-hadron decays. First, we explicitly veto contributions from B~ —» D°D~, with 
either D~ —» 7r - 7r + 7r - or D~ —y K~7T + n~, by rejecting candidates in which the X~ system 
has invariant mass within 15MeV/c 2 of the known D~ mass. Contamination from other 
final states that include a charmed particle are also sought by forming all two-, three- and 
four-body combinations (except the D —y h + h'~ signal decay), and checking for peaks at 
any of the known charmed particle masses. Contributions from D° —> K~n + , K~tt + 71~tt + , 
Df —y K + K~n + and D + —y K~n + 7i + decays are seen, and ±15MeV/c 2 mass vetoes are 
applied around the known charm particle masses. In addition, D* + contributions are 
removed by requiring the invariant mass difference, M[{K~ti + )£,ti + ] — M[(K~tt + )d\ > 
148.5 MeV/c 2 . This removes both partially reconstructed B —>• D* + X final states and fully 
reconstructed states, such as B —» D 1 (2420)°/r - , D 1 (2420)° —* D* + n~, D* + —> D°n + 
signal decays. The latter, while forming a good signal candidate, are flavor-specific, 
and therefore would reduce the coherence of the final state. Those D**° —y D°7r + 7r _ 
contributions that do not have a D* + intermediate state are kept, since they are not 
flavor- specific. 

Another potentially large source of background is from five-body charmless B decays. 
Unfortunately, their branching fractions are generally unknown, but they are likely to be 
sizable compared to those of the B~ —$■ DX~ signal decays. Moreover, these backgrounds 
could have large CP asymmetries, as seen in three-body 5-meson decays 121,39,40 . It is 
therefore important to suppress their contribution to a negligible level. This is investigated 
by applying all of the above selections, except that D candidates are selected from a D mass 
sideband region instead of the signal region. The sideband region is chosen to avoid the 
contribution from the other two-body D decays with one misidentihed daughter. Charmless 
backgrounds are seen in all modes. These backgrounds are reduced to a negligible level by 
requiring that the D decay vertex is displaced significantly downstream of the B~ decay 
vertex, corresponding to three times the uncertainty on the measured D decay length. A 
more stringent requirement, corresponding to five times the uncertainty on the measured 
D decay length, is imposed on the B~ —» [n + Ti~] D Xj d decays, which is found to have a 
much larger charmless contribution. After these requirements are applied, the charmless 
backgrounds are consistent with zero, and the residual contribution is considered as a 
source of systematic uncertainty. 

Another important background to suppress is the cross-feed from the ADS CF B~ —» 
[K~n + ] D X- decay into the ADS DCS B~ —» [K + tt^]bX~ sample, which may happen 
if the K~ and n + are both misidentihed. Since the CF yield is expected to be several 
hundred times larger than that of the DCS mode (depending on the values of r#, 5b, k 
and 7 ), a large suppression is necessary. The combined D° mass and PID requirements 
provide a suppression factor of 6 x 10 -5 . An additional requirement that the Kir invariant 
mass (after interchanging the K~ and n + masses) differs by at least 15 MeV/c 2 from the 
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known D° mass decreases the suppression level to 0.9 x 10 -5 . This leads to a negligible 
contamination from the CF ADS mode into the DCS decay. The same veto is applied 
to both the ADS CF D° —* K~n + and DCS D° —y K + ir~ decays, so that no efficiency 
correction is needed for R x . 

Lastly, in order to have a robust estimate of the trigger efficiency for signal events, we 
impose requirements on information from the hardware trigger; either (i) one or more of the 
decay products of the signal candidate met the trigger requirements from the calorimeter 
system, or (ii) the event passed at least one of the hardware triggers, and would have 
done so even if the signal decay was removed from the event. These two classes of events 
constitute about 60% and 40% of the signal candidates, respectively, where the overlap is 
assigned to category (i). 

The selection efficiencies as a function of several two- and three-body masses in the 
B~ —y [K~7i + ]£)Xj decay are shown in Fig. [lj for both the GLW and ADS selections. 
The efficiencies for other D final states are consistent with those for D —» K~tt + . The 
m(Dir-) and m(7r + 7r - ) efficiencies include two entries per signal decay, as there are two 
7 t~ in the final state. The analogous efficiencies for the B —>■ DX~ decay are shown in 
Fig-! The relative efficiencies of the ADS to GLW selections are consistent with being flat 
across each of these masses. These efficiencies include all selection requirements, including 
PID. However, events in which any of the signal decay products is outside of the LHCb 
detector acceptance are not included, since they are not simulated; thus to obtain the total 
selection efficiency, these efficiencies should be scaled by a factor of 0.11, as determined 
from simulation. 

Figure [ 3 ] shows the A% and X~ invariant mass distributions for B~ —> [K~tt + }oXj 
and B~ —> [K~tt + ]dX~ signal decays after all selections, except for the X~ and K*° 
mass requirements. These signal spectra are background subtracted using the sPlot 
method |4l], with the B candidate invariant mass as the discriminating variable. The X^ 
and X~ contributions peak in the region below 2 GeV/c 2 , consistent with the dominance 
of resonances such as ai(1260)~ —> 7r _ 7r + 7r _ to the Xj system, and one or more excited 
strange resonances contributing to X~. The dip at 1.97 GeV/c 2 is due to the Dj mass 
veto. 


5 Fits to data 

The signal yields are determined through a simultaneous unbinned extended maximum 
likelihood fit to the 16 B ± candidate invariant mass spectra. These 16 spectra include the 
four B~ —y DX^ decays, where D —y iF ± 7T =F , K + K~ and 7r + 7r _ , the corresponding four 
charge-conjugate decays, and the set of eight modes where Xj is replaced with X~. The 
signal and background contributions across these modes are similar, although not identical. 
Where possible, common signal and background shapes are used; otherwise simulation 
is used to relate parameters in the lower yield modes to the values obtained from the 
high yield CF D —> Kn modes. Signal and background yields are all independent of one 
another in the B + and B~ mass fits; thus CP violation is allowed for all contributions in 


7 



m(7t + 7t') [MeV/c 2 ] 


m(Dr) [MeV/c 2 ; 


m(D7t + ) [MeV/c 2 ] 




+GLW 

+ADS 


Figure 1: Signal efficiencies for the B —»• [K 'K + ]£ ) X d decay when applying the GLW and ADS 
selections. The efficiencies are shown as a function of five different two- and three-body masses. 


the mass spectrum. Unless otherwise noted, the shapes discussed below are obtained from 
simulated decays. 


5.1 Signal shapes 

The B~ mass signal shapes are each parameterized as the sum of a Crystal Ball (CB) 
shape 42 and a Gaussian (G) function, 


F,ig oc /cBCB(m B ,cr C B,acB,n) + (1 - /c B )G(m s , <r„). 


(13) 


The Gaussian function accounts for the core of the mass distribution, whereas the CB 
function accounts for the non-Gaussian radiative tail below, and a wider Gaussian resolution 
component above, the signal peak. A small difference is seen between the shapes for the 
B~ — > DX d and B~ —> DX~ decays, and so a different set of signal shape parameters 
is used to describe each, except for a common value of the fitted B~ mass, ms- The 
signal shapes are not very sensitive to the power-law exponent, n, which is fixed to 
10. The parameters acB; & g an d /cb are allowed to vary freely in the fit. to the data. 
From simulation, we find that for all 16 modes, Ocb/cG consistent with 1.90, and 
this ratio is imposed in the fit. Simulation is also used to relate the mass resolution in 
the D — > Ii + K~, 7 r + 7 modes to that of the D — >■ K it mode, from which it is found 


8 











m(7t + 7r) [MeV/c 2 ] m(KV) [MeV/c 2 ] m(Dr) [MeV/c 2 ] 



+GLW 

+ADS 


Figure 2: Signal efficiencies for the B —»• [K -k + ]£)X s decay when applying the GLW and ADS 
selections. The efficiencies are shown as a function of five different two- and three-body masses. 


that af K]oX = (0.947 ±0.011)4^ 1dA ' and <t [ ™ ] dX = (1.043 ± 0.011)(t^ ] dA .The 
relations are consistent between the B~ —y DX d and B~ —> DX~ modes, and are applied 
as fixed constraints (without uncertainties) in the mass fit. 

5.2 Backgrounds and their modeling 

The primary sources of background in the mass spectra are partially reconstructed B —> 
D&X~ decays, cross-feed between B —> DX d and B —» DX S , and other combinatorial 
backgrounds. All of the spectra have a contribution from combinatorial background, the 
shape of which is described by an exponential function. Its slope is taken to be the same 
for the CP-conjugate B~ and B + decays, but differs among the various D and X~ final 
states. 

The main contribution to the partially reconstructed background comes from B~ —» 
[P°7r°, D° r y] D *oX~ or B° —> [D 0, n + ]^*+ X~ decays, where a pion or photon is not considered 
when reconstructing the B~ candidate. Because the missed pion or photon generally has 
low momentum, these decays pass the full selection with high efficiency. The shapes of 
these distributions are modeled using parameterized shapes based on simulated decays. 
Since the Dalitz structure of these backgrounds is not known, we do not rely entirely on 
simulation to reproduce the shape of this low-mass component. Instead, the parameters 
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Figure 3: Signal distributions of the (left) XT invariant mass in B~ —> DXT decays and (right) 
X~ invariant mass, in B~ —> DX~ decays, for D —> K~n + . The distributions are obtained 
using the sPlot method. In both cases, all selections, except the M{X~) < 2GeV/c 2 and the 
K*° mass selection, are applied. The dip at 1.97 GeV/c 2 is due to the Df meson veto. 


of the shape function that depend on the decay dynamics are allowed to vary freely, 
and are determined in the fit. The shape parameters for these backgrounds are varied 
independently for B~ —>■ DX^ and B~ —>• DX~ decays. 

Another background contribution which primarily contributes to the B~ —> DX~[ ADS 
suppressed mode is the B° —> D°7 r _ 7r + 7r _ 7r + decay, where there is no D* + intermediate 
state. This decay can contribute to the ADS CF mode if a tt + is excluded from the 
decay, or to the ADS DCS mode if a 7r~ is not considered. The branching fraction 
for this decay is not known, but the similar CF decay B° —> D*°7r~n + 'n'~7r + is known 
to have a relatively large branching fraction of (2.7 ± 0.5) x 10~ 3 43,44 . Assuming 


B(B° -> D°ir~ 


7T 7T 7T 


~ B(B° -)• 


D*°7T 7T + 7T 


7T 1 


this background contribution is 


about two orders of magnitude larger than the DCS signal, although it peaks at lower 
mass than the signal. The selection efficiency and shape of this background are difficult 
to determine from simulation, since there have not been any studies of this final state 
to date. Its shape is obtained from simulations that assume a quasi two-body process, 
B° —y D°R, R —y 7r _ 7r + 7r _ 7r + , which decays uniformly in the phase space. An ARGUS 
shape 45 convolved with a Gaussian function provides a good description of this simulated 
background. Its shape parameters are shared between B + and B~ and are allowed to 
vary freely in the fit, except for the Gaussian width, which is fixed to the expected mass 
resolution of 15 MeV/c 2 . 

The analogous B° —> D 0 K~7T + n~'n + decay does not pose the same contamina¬ 
tion to the DCS ADS B + —» [K~tt + ] d X+ signal, since a missed 7r~ leads to a 
B + —y D°K~7 r + 7r + candidate, which is not one of the decays of interest. However, 
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in the B Q S —> [K ti + ]d° 


K + 7T 


TT^TT 


presence of the s quark within the meson. This decay is unobserved, but the similar 


decay, opposite-sign kaons are natural due to the 
o 

_ S 

decay, B° s —> D 0 K + n~ 1 has a relatively large branching fraction of (1.00 ±0.15) x 10~ 3 46 
Based on other 5-meson decays, one would expect the B° s —> D°K + 7 t~7t + tt~ decay to be 
at the same level, 0(1O -3 ), which is two orders of magnitude larger than the signal. The 
shape of this background has a similar threshold behavior as for the B° —> 
decay discussed previously, and therefore its contribution is also modeled from simulated 
decays using an ARGUS shape convolved with a Gaussian function with freely varying 
shape parameters. 

In the fit, we also model cross-feed between the B~ —>■ D *' 1 Xj and 5~ —>■ DX~ 
decays. The shapes of these cross-feed backgrounds are obtained from simulation. The 
cross-feed rate is obtained from D* + —)■ D°tt + , D° —> K~n + calibration data, reweighted 
to match the properties of the signal decays. All selection requirements on the B~ —>■ DX~ 
decays, including |M(A' _ 7t + ) — M K * o| < 100MeV/c 2 and M(X~) < 2 GeV/c 2 , are taken 
into account. In total, we find that 0.66% of 5~ —>■ DX^ are misidentified as B~ —s- DX~ 
for the GLW modes and 0.16% for the ADS modes. The lower value for the ADS modes 
is due to the tighter PID requirements on the K~ candidate in the X~ system. The 
cross-feed from B~ —s- DX~ into B~ —s- DX^ is evaluated in an analogous manner, and 
is found to be 13.7%. Since the ratio of branching fractions is B(B~ —> DX~)/B(B~ —$■ 
DXj) ~ 0.09 |47|, the yield of this background is only about 1% of the signal yield. 

Other sources of background that contribute to the B~ —y DX~ modes are the 
B~ — > D°[K~K + tt~] d - and B~ — > D°K~[K + tt~]^* 0 decays, where the K + is misidentified 
as a 7r + meson. The shapes are similar for these two backgrounds and thus a single shape 
is used, based on a parameterization of the B~ candidate mass distribution in simulated 
B~ —>• D°[K~K + n-] D - decays. Taking into account known branching fractions 21 

->■ dV 


efficiencies from simulation, and K + —» n + misidentification rates from D* 
calibration data, we expect a contribution of 1.6% of the B~ —» DX~ signal. 


5.3 Fit results 

The invariant mass spectra for the 5~ —> DX~ ADS and GLW signal modes are shown in 
Figs. [4] and [ 5 J with the corresponding spectra for the B~ —y DX~[ normalization modes 
in Figs. [6] and [TJ Results from the fits are superimposed along with the various signal 
and background components. The fitted yields in the ADS and GLW modes are given in 
Tables |T] and [2] 

Highly significant signals are seen in all modes, except for the ADS DCS B~ —> DX~ 
decay. This is the first time these decays have been observed in modes other than the CF 
D° —>■ K~tt + decay. Figure [8] shows the suppressed ADS mode, B ± —>• D[K + 7r^} D K ± 7r lz 7r ± : 
summed over both 5-meson charge states. The significance of the peak, which exceeds 
three standard deviations, is discussed later. 
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Table 1: Fitted yields in the ADS modes with / = Kir, for the signal and corresponding 
normalization modes. 


Decay mode 

B yield 

B + yield 


(N f ) 

AW 

B ± ->■ DXf, D -)• K~tt + 

36 956 ± 214 

37843 ±219 

B ± -)• DXf, D -► K + 7f- 

161 ± 20 

162 ± 20 


(N f ) 
y at,xr J 

AW 

B ± -)• DXf, D K~ti + 

1234 ± 37 

1226 ± 37 

B ± -)• DXf, D -> K + n~ 

13.0 ±5.3 

6.6 ±4.0 


Table 2: Fitted yields used in the GLW analysis with / = K±ir ±, K + K and 7r + 7r , for the 
signal and corresponding normalization modes. 


Decay mode 

yield 

AW 

5+ yield 

AW 

B ± -)• DXf, D -> K~ 7T+ 
W -)• D -> K + K~ 

B ± ->• DA±, D 7T+7T- 

45 213 ±226 

3899 ± 63 

1669 ± 38 

46 488 ± 230 
4084 ± 65 
1739 ± 40 


(A 7 ) 

AW 

W -)• D -► A-vr + 

W -)• DAf, D -)■ K + K~ 
B ± DAf, D ->• TT+Tr" 

1699 ±47 
155 ± 14 

59 ±9 

1744 ± 47 
171 ± 14 

70 ±9 


6 Determination of CP observables 


The CP observables are obtained by expressing the fitted signal yields in terms of corrected 
yields and the CP parameters. For the decay B ± —> fpXf, where fp is either the ADS 
CF decay or a CP eigenstate, the fitted yields can be written as 


n! = - 


N, 


f 


corr,Xd 


fit,x2 


i _i_ pf 

1 + v*. 


( x ± ^Lw,x d ) + C j. Xi 


±5 

d 


(14) 


where Nf OTrXd is the total corrected yield (sum of B and B + ), F^ x are the estimated 


fractions of signal events removed by the D° and Z)j*j + vetoes, (±f Y± are the estimated 


charmless background yields, and Aj., dW Xd is the raw CP asymmetry. 

The htted yields in the corresponding B ± —> DXf decays are written in terms of 
the corrected B ± —> DXf yields in Eq. 14 and the CP observable R[/ d defined in Eqs. |9j 
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Figure 4: Mass distributions of B~ —» DX~ candidates using the ADS selections, for (top left) 
B~ —f [K~tt + ]dX~ , (top right) B + —> [K + tt~]dX^, (bottom left) B~ —> [K + tt~]£)X~ , and 
(bottom right) B + —> [K~tt + ]dX+ . 


and 10 


as 


nL 


N U = §<< {) (1 T aL - xJ + c l*f 


(15) 


where the meaning of the symbols parallels those in Eq. 14 


For the ADS suppressed modes, the four DCS yields are expressed in terms of 
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Figure 5: Mass distributions of B~ — > DX~ candidates using the GLW selections, for (top 
left) B~ -* [K~tt + ]dX~ , (top right) B + —> [K + tt~}dX+, (middle left) B~ —> [K + K~]dX~, 
(middle right) B + —> [K + K~]dX+ , (bottom left) B~ —> [7r + 7r _ ]£)X“, and (bottom right) 
B + [7r + 7r _ ] D X+. 
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Figure 6: Mass distributions of B — 
B~ [K~tt + ] d Xj, (top right) B + 
(bottom right) B + —> [K~tt + ]dX^. 


d candidates using the ADS selections, for (top left) 


DX u 

-> [K+7T-] D Xt 


d ’ 


(bottom left) B —y [ I\ + tt ]nX, , and 


the corrected CF yields, N^ rr 


corr ,Xd 5 


as 


N^ 1 1 

&t,X± 


N KTw ± 


- I 

— I traw 




= IB. 


N, 




corv,Xd 


1 + f $£ 


N. 


K± 71^ 


corr,X s 


1 _J_ pK Ttt± 

1 + Vas 


I ,^yA'T 7T 1 * 1 

+ VJ . 


+ c t,xt • 


( 16 ) 

(U) 
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Figure 7: Mass distributions of B~ —>• 17X)" candidates using the GLW selections, for (top 
left) B~ —> [K~tt + \dX^, (top right) B + —> [I\ + tt~]dX ^, (middle left) B~ -» [A' + Ar~]£>X)", 
and (middle right) H + -» (bottom left) F> - -» [7r + 7r _ ]i)A 7 ' ( )", and (bottom right) 

B+-> [tt+it-} d X+ 
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Figure 8: Mass distributions for the suppressed ADS mode, B± 
B + and B~). 


(sum of 


where (1 =)= Af^x s ) gives the corrected yield for the favored B ± —> 

[K ± 'K T ] D xf decays. 

The corrections for the D° and D { '*\ + vetoes, FL „ , are determined by interpolating 
from the mass regions just above and below the veto region, and lead to corrections that 
range from 0.6% to 5.8% of the expected yield. Uncertainties on these corrections are 
considered as sources of systematic uncertainty. Potential contamination from charmless 
five-body decays is determined by fitting for a B ± signal component when the D candidates 
are taken from the D° mass sideband region, as described previously. The charmless 
contributions are negligible, and the uncertainties are included in the systematic error. 


The yields, as determined from the fitted values of the CP parameters in Eqs. [14]fT7 
given in Tables [I] and [2j 


are 


The raw observables, *4.(. aw Y and R^ aw include small biases due to the production 
asymmetry of B ± mesons, Ab± (affecting *4.f aw x only), and from the detection asymmetries 
of kaons and pions, Ak and A n . The corrected quantities are then computed according to 
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a k+k ~ 

^x d 

= ^-raw^Y d — Ab± — A w , 

(18) 

A 7T + 7T — 

^X d 

= AawPQ — Ab± — A n , 

(19) 

aK~ 7T+ 

^x d 

= ALvPQ ~ Ab± — Ak , 

(20) 

A K+K ~ 

x 3 

= — Ak , 

(21) 

XxX 

= Aawpq, — Ab± — Ak 1 

(22) 

A K~ 7T+ 

= -^rawW s — Ab± — 2Ak + Arj 

(23) 

R x t 

= R*i( 1 - 2 Ak + 2A), 

(24) 

R x * 

= R*i(l + 2A K - 2A n ), 

(25) 

R x » 

= r£( 1 - 2 Ak + 2A), 

(26) 

R x > 

= R x i(l + 2A K -2A n ). 

(27) 


The pion detection asymmetry of An = 0.000 ± 0.003 is obtained by reweighting the 
measured A detection efficiencies 48 with the expected momentum spectrum for signal 
pions. The kaon detection efficiency of Ak = —0.011 ± 0.004 is obtained by reweighting 
the measured K~n detection asymmetry (49i| using the momentum spectrum of signal 
kaons, and then subtracting the above pion detection asymmetry. For the production 
asymmetry, the value Ab± = —0.008 ± 0.007 is used 50 , based on the measured raw 


asymmetry in B ± —> J/'^>K ± decays 51 and on simulation. 


6.1 Systematic uncertainties 

Most potential systematic uncertainties on the observables are expected to cancel in 
either the asymmetries or ratios that are measured. The systematic uncertainties that 
do not cancel completely are summarized in Table [3] The PID and trigger asymmetries 
are evaluated using measured kaon and pion efficiencies from D* + —> D°7t + calibration 
samples in data that are identified using only the kinematics of the decay. The efficiencies 
for the B + and B~ signal decays are then obtained by reweighting the kaon and pion 
efficiencies using simulated B ± —> DX ± decays to represent the properties of signal data. 
We find no significant charge asymmetry with respect to the PID requirements, and use 
= 0.000 ± 0.006, where the uncertainty is dominated by the finite sample sizes of 
the simulated signal decays in the reweighting. The asymmetry of the hardware trigger is 
assessed using measured hadron trigger efficiencies in D* + —> D°ir + , D° —> K~tt + decays, 
reweighted to match the momentum spectrum of tracks from signal decays. Defining the B ± 
hadron trigger efficiency as €b±, the charge asymmetry of the trigger (e#- — e_B+)/(e_B- +^s+) 
varies from 0.000 ± 0.003 for B~ —>■ [K + K~] D X~ to 0.007 ± 0.003 for B~ —> [K + n + ] D X~. 
These values are applied as corrections. 

On Rcp+ and R x , we have either a double-ratio or a ratio of final states with 
identical particles (apart from the charges), and therefore there is a high degree of 
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Table 3: Systematic uncertainties, in percent, on the fitted parameters. 


Source 

D 

A(B ± DXf) 
h + h~ Kit 

A(B ± DXf) 
h + h~ Kn 

Rcp+ 

K+K~ 7T+7T- 

R* 

Kn 

^ a 

Ab± 

0.7 

0.7 

0.7 

0.7 

- 

- 

- 

- 

Ak 

- 

0.4 

0.4 

0.8 

- 

- 

0.7 

0.7 

An 

0.3 

- 

- 

0.3 

- 

- 

0.6 

0.3 

Trigger 

0.4 

0.4 

0.4 

0.4 

1.5 

1.5 

1.5 

1.5 

PID 

0.6 

0.6 

0.6 

0.6 

1.2 

1.2 

1.2 

1.2 

Signal model 

- 

- 

- 

- 

1.1 

1.1 

- 

- 

Bkgd. model 

- 

- 

- 

- 

1.6 

1.6 

4.0 

10.0 

Charmless back. 

- 

- 

- 

- 

1.0 

1.0 

1.0 

1.0 

Cross-feed 

- 

- 

- 

- 

1.0 

1.0 

1.0 

1.0 

D vetoes 

- 

- 

- 

- 

1.0 

1.7 

1.0 

1.0 

Rcp+ approx. 

- 

- 

- 

- 

1.0 

1.0 

- 

- 

Total 

1.0 

1.1 

1.1 

1.3 

3.4 

3.8 

4.9 

10.4 


cancellation of potential systematic uncertainties. We expect that for these ratios, the 
relative trigger efficiencies would yield a value close to unity. After reweighting the 
measured trigger efficiencies according to the kinematical properties of signal decays 
(obtained from simulation), we find that the ratios of trigger efficiencies are within 1.5% 
of unity, which is assigned as a systematic uncertainty. Using an analogous weighting 
procedure to the measured PID efficiencies, we find that the relative PID efficiency is 
equal to unity to within 1.2%, which is assigned as a systematic uncertainty. 

We also consider uncertainty from the signal model, the background model, the charm¬ 
less contamination, the D vetoes, and the detection asymmetries. For the signal model 
uncertainty, all of the fixed signal shape parameters are varied by one standard deviation, 
and the resulting changes in the CP parameters are added in quadrature to obtain the total 
signal shape uncertainty (1.1%). For the background-related uncertainties, we consider 
a polynomial function for the combinatorial background, and vary the fixed background 
shape parameters of the specific 6-hadron backgrounds within their uncertainties, and add 
the deviations from the nominal result in quadrature (1.6%). For the ADS-suppressed 
modes, larger uncertainties are assigned based on an incomplete understanding of the 
contributions to the low mass Bf^ —> D°X background. 

The charmless backgrounds are all consistent with zero, and the uncertainty is taken 
from fits to the D sideband regions (1.0%). Uncertainties due to the cross-feed contribu¬ 
tions (such as B~ —> DXj reconstructed as B~ —> DX ~) are assessed using simulated 
experiments, by simulating the mass distributions with a larger cross-feed and fitting with 
the nominal value (1.0%). The uncertainties due to vetoing potential contributions from 
other D mesons are assessed by interpolating the mass spectrum just above and below the 
veto region into the veto region. The associated uncertainties are all at the 1.0% level, 
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except for the the B —» [7r + 7r ]dX s mode, which has an uncertainty of 1.7%. 

The uncertainties on the ratios Rgjj 1 ' and R Xs ’ d are each summed in quadrature, 
giving total uncertainties in the range of (3.4 — 10.4)%, depending on the mode. 


7 Results and summary 


The resulting values for the CP observables are 


Rcp+~ = 1-043 ± 0.069 ± 0.034, 

R*cpX = 1-035 ± 0.108 ± 0.038, 

A x+ d K ~ = -0.019 ± 0.011 ± 0.010, 

Ax'*' = -0.013 ± 0.016 ± 0.010, 

A x ~ n+ = -0.002 ± 0.003 ± 0.011, 

R x d = (43.2 ± 5.3 ± 2.1) x 10~ 4 , 

R x d = (42.1 ± 5.3 ± 2.1) x 10~ 4 , 

A x+ s K ~ = -0.045 ± 0.064 ± 0.011, 

Ax'*' = -0.054 ± 0.101 ± 0.011, 

Ax~* + = 0.013 ± 0.019 ± 0.013, 

R xt = (107^44 ± 11) x 10~ 4 [ < 0.018 at 95% CL ], 
R X7 = (53 tf 2 ± 6) x 10 -4 [ < 0.012 at 95% CL ]. 

The values of Rcp+ are averaged to obtain 


R CP+ = 1.040 ± 0.064, 

where the uncertainty includes both statistical and systematic sources, as well as the 
correlations between the latter. 

The significances of the suppressed ADS modes are determined by computing the ratio 
of log-likelihoods, \J2 log(£ 0 /£ min ), after convolving C with the systematic uncertainty. 
From the value of C at the minimum (£ m j n ), and the value at R Xs = 0 (/%), the 
significances of the non-zero values for R Xs and R Xs are found to be 2.On and 3.2<r, 
respectively. The overall significance of the observation of the ADS suppressed mode is 
obtained by adding the log-likelihoods, resulting in a significance of 3.6 standard deviations. 
This constitutes the first evidence of the ADS suppressed mode in B~ —> DK~tt + tt~ 
decays. 

For completeness, we also compute the related observables Rads an d ^4ads, which are 
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commonly used. For the B —y DX S modes, the values are 


R 


x s _ 

ADS 


= (R x ° + R x ° )/2 = ( 85 I 33 ) x 1(T 4 , 




+36 \ 


R 


x: 


R 


•xi 


A A » = 

ADS R x * + R x ° 


_ _n OO+0.36 

— U.OO_ 0 34. 


For the favored modes, the corresponding values are 


i?ADS = ( R + R )/ 2 = (42.7 ± 5-6) x 10 


-4 


A Xd = 

^ADS — 




R x ° -R 

rx- + R x; 


= -0.013 ±0.087. 


The averages are computed using the asymmetric uncertainty distributions, and include 
both statistical and systematical sources. 

To assess the constraints on 7 that these observables provide, they have been imple¬ 
mented in the litter for 7 described in Ref. 


14 . Two fits are performed, one that uses 


only information from B —> 
B~ —y DX~ and B —y DX, 


DX s , and a second that uses the observables from both 


system, ro, hp 71 


decays. In both fits, the parameters from the D-nreson 

are constrained in an analogous 

and B~ 


Af P (K+K ~), and Afl r P 


-+ DK 


yr 1 7T 

—> Dn~ 


case 


14 . The four parameters 


%d, Ud, 

way to what was done for the B 
rp, 8 b, k and 7 are freely varied in each fit. In the combined fit, three additional strong 
parameters, r p Xd , 8 p Xd , n DXd are included, which are analogous to those that apply to 
the B~ —* DX~ decay. 

The projections of the fit results for 7 , rs and tb versus 7 , are shown in Fig. [9] using 

The value of 7 is found to be ( 74 ^ 3 )° for the 
—» DX~ and B~ —> DX 


52 


14 


(see also Refs 

only fit, and ( 74 ^ 49 )° for the for the combined B 


d 


the method of Ref. 

B~ DX, 

fit. The value of Tb is nearly identical in the two cases, with corresponding values of 
r B = 0 . 081 !q'q 27 and r B = O.OSllg^g- ^- s expected, most of the sensitivity comes from 
the B~ —>• DX~ decay mode. This value is almost identical to the LHCb combined result 
of ( 73 ^ 4 o)° found in Re: 

B~ —> DK~ decays 1:50 
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56 


The value of re is similar to the values found in other 

found in 


but smaller than the value of 0.240^ o 48 
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neutral R-meson decays. The strong phase 8b, averaged over the phase space, peaks at 
172° for both fits, but at 95% CL all angles are allowed. The constraints on the coherence 
factor are relatively weak; while the most likely value is close to 1 , any value in the interval 
[ 0 , 1 ] is allowed at one standard deviation. 

In summary, a pp collision data sample, corresponding to an integrated luminosity of 
3.0 fb -1 , has been used to study the B~ —y DX~ and B~ —y DX^ decay modes, where the 
D meson decays to either the quasi-flavor-specific Kn final state or the K + K~ and 7 r + 7 r _ 
CP eigenstates. We observe for the first time highly significant signals in the CP modes 
for both the favored and suppressed B~ decays, and we also report the first evidence for 
the ADS DCS B~ —» [K + 7t^]dK~tt + tt~ decay. We measure the corresponding ADS and 
GLW observables for the first time in these modes. A fit for 7 using only these modes 
is performed, from which we find 7 = (74^23)° for the fit with only B~ —> DX~, and 
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Figure 9: Projections of 1—CL versus (left) 7, (right) r B ~+ DK ^ ^ 5 and (bottom) 
r B ~> DK 7r+7r versus 7, using only B~ — > DK~tt + tt~ decays, and the combination of 
B~ — y DK~tt + tt~ and B~ —> i>7r~7r + 7r _ decays. The 68.3% and 95.5% confidence level 
(CL) limits are indicated for the 7 and rs projections. The 39% level contours in r B 
versus 7 correspond to the 68.3% level contours in the one-dimensional projections. 




7 = (74+^)° for the combined B~ —> DX~ and B~ —> DX^ fit. Values of 7 below about 
25° and larger than approximately 165° are not excluded by these modes alone, but are 
excluded when other modes are considered 114]. The precision on 7 in this analysis is 
comparable to, or better than, most previous measurements. 
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